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Abstract 

Recent developments of high-order CCM have been to extend existing formalism and codes to 
s > \ for both the ground and excited states. In this article, we describe how "generalised" expec- 
tation values for a wide range of one- and two-body spin operators may also be determined using 
existing the CCM code for the ground state. We present new results for the spin-spin correlation 
functions of the spin-half square- and triangular-lattice antiferromagnets by using the LSUBm ap- 
proximation. We show that the absolute values of the spin-spin correlation functions |(s(0).s(r))| 
converge with increasing approximation level for both lattices. We believe that the LSUBm ap- 
proximation provides reasonable results for the correlation functions for lattice separations roughly 
of order r ~ m for the square lattice. We compare qualitatively our results for the square lattice 
to those results of quantum Monte Carlo (QMC) and we see that good correspondence is observed. 
Indeed as seen by QMC, the spin-spin correlation function initially decays strongly with |r| before 
becoming constant for larger values of \r\. CCM results are also compared to results of exact 
diagonalisations for both lattices. ED results demonstrate a strong finite-size effects at lattice sep- 
arations r = L/2 (where N = L x L) for both lattices. The behaviour of the correlation function 
for the triangular lattice is qualitatively similar to that of the square lattice, namely, that it decays 
strongly at first before becoming constant. This is in keeping with the behaviour of both models, 
which are believed strongly to be Neel-ordered from approximate studies. The CCM has been 
shown many times to provide consistently good results for the ground-state energy and sublattice 
magnetisation of a wide range of quantum spin models. Here we have shown that the CCM also 
provides good results for the spin-spin correlation function. 
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I. INTRODUCTION 



The coupled cluster method (CCM) is a well-known method 

of quantum many-body theory (QMBT). The CCM has been applied with much success 



in order to stu dy quan tum magnetic 
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52j). In particular, the use 
52| | of the CCM has been found to be 



very effective with respect to these spin-lattice problems. Recent developments of high-order 
CCM formalism and codes have been to treat systems with spin quantum number of s > ~ 
for both the ground and excited states 52| . In this article, we show how these ground-state 
formalism and codes may also be used directly to find "generalised" expectation values; that 
is, expectation values for a wide range of one- or two-body spin operator that are defined 
prior to the CCM calculation. 

We apply the CCM to the spin-half Heisenberg model on the square and triangular 
lattices. The Hamiltonian is specified as follows, 



H 



"3 > 



(1) 



where the sum on (i, j) counts all nearest-neighbour pairs once. The ground states of all of 
the cases considered here are classically ordered, albeit by a reduced amount due to quantum 
fluctuations. Indeed, the best estimates of the amount of classical ordering of the square 
lattice case from approximate methods are 61% from CCM 4l|], 61.4% from quantum Monte 
Carlo studies 531] . 61.4% from series expansions 54j . 61.38% from spin- wave theory 551] . and 



63.4% from exact diagonalisations 56j] . Good results for the spin-spin correlation function 



of the spin-half square-lattice Heisenberg antiferromagnet were found using quantum Monte 
Carlo in Ref. 5tJ on lattices of size up to N = L x L, where L = {16, 22, 32}. They observed 
that the correlation functions decayed with separation r. However, they also saw finite-size 
effects, indicated by cusp-like behaviour in these correlation functions, at distances given by 
L/2. Exact diagonalisations have also been carried out for the correlation function of the 
square lattice 58] . For the spin-half triangular-lattice antiferromagnet, the best estimates of 
the amount of classical ordering are 41% from quantum Monte Carlo studies 59j, 40% from 



series expansions 



601 ] . 39% from exact diagonalisations 56j , and 40% from previous CCM 
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calculations 41|. Very few results for the spin-spin correlation functions of the triangular 
lattice seem to exist apart from results of exact diagonalisations (e.g., see 0]). 



II. METHOD 



The details of the practical application of high-order coupled c 



uster method (CCM) 



formalism to lattice quantum spin systems are given in Refs. [22j, [26|, [32|, [38|, [52| and also in 
the appendices to this article. However, we point out now that the ket and bra ground-state 
energy eigenvectors, and of a general many-body system described by a Hamiltonian 
H, are given by 

H\9)=E g \9); (9\H = E g {9\. (2) 

Furthermore, the ket and bra states are parametrized within the single-reference CCM as 
follows: 



(§| = (<S>\Se~ s ; S = 1 + J2 Si C 7 ■ 



(3) 



One of the most important features of the CCM is that one uses a single model or reference 
state |$) that is normalized. We note that the parametrisation of the ground state has the 
normalization condition for the ground-state bra and ket wave functions = ($1$) = 

1). The model state is required to have the property of being a cyclic vector with respect to 
two well-defined Abelian subalgebras of multi-configurational creation operators {Cf } and 
their Hermitian-adjoint destruction counterparts {Cj = (C/ )'}. The interested reader is 



521 ] for more information regarding how the CCM 



referred to the Appendices and to Ref. 
problem is solved for. 

Here, we use the classical ground state as the model state. For the square lattice, this is 
the Neel state in which neighbouring spins are antiparallel, and, for the triangular lattice, 
this a Neel-lke state in which neighbouring spins on three sublattices are at 120° to each 
other. For the square lattice, we perform a rotation of the local axes of the up-pointing spins 
by 180° about the ?/-axis. The transformation is described by, 



s\ s z 



(4) 
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The model state now appears mathematically to consist of purely down-pointing spins. In 
terms of the spin raising and lowering operators = s% ± is v k the Hamiltonian may be 
written in these local axes as, 

l N r i 

^ = "oE sis+ + srsj + 2s z i s z j , (5) 

where the sum on again counts all nearest-neighbour pairs once on the square lattice. 
Again, the classical ground-state of the Heisenberg model of Eq. (ED) for the triangular 
lattice is the Neel-like state where all spins on each sublattice are separately aligned (all in 
the xz-plane, say). The spins on sublattice A are oriented along the negative z-axis, and 
spins on sublattices B and C are oriented at +120° and —120°, respectively, with respect to 
the spins on sublattice A. We again rotate the local spin axes of those spins on the different 
sublattices. Specifically, we leave the spin axes on sublattice A unchanged, and we rotate 
about the y-axis the spin axes on sublattices B and C by —120° and +120° respectively, 

x 1 x z x 1 x a/3 z 

s b ~> ~2 Sb 2~ Sb ' S ° ~~ ^ ~2 S ° ~2~ Sc ' 

Sy v c*y ■ c*y ^ 
B *S j b C ~* S C i 

z v^3 x 1 z z V^3 x 1 z , . 



Once again, the model state now appears to consist of purely down-pointing spins. We may 
rewrite Eq. (TjQ) in terms of spins defined in these local quantisation axes for the triangular 
lattice, such that 



(i->i> 



1 g g V^^" / ~ I Z — I ~ _ 2, 

~2~ S i S j + ~ \ S i S j + S i S j ~ S i S j ~ S i S j 



1 3 

+s {4*7 + »t4) ~ § (44 + s i4) } • ( ? ) 

We note that the summation in Eq. ([7]) again runs over nearest-neighbour bonds, but now 
also with a directionality indicated by (i — > j), which goes from A to B, B to C, and C to A. 

The CCM formalism is exact in the limit of inclusion of all possible multi-spin cluster cor- 
relations within S and S, although this is usually impossible to achieve practically. Hence, 
we generally make approximations in both 5" and S. The three most commonly employed 
approximation schemes previously utilised have been: (1) the SUBn scheme, in which all 
correlations involving only n or fewer spins are retained, but no further restriction is made 
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concerning their spatial separation on the lattice; (2) the SUBn-m sub-approximation, in 
which all SUBn correlations spanning a range of no more than m adjacent lattice sites are 
retained; and (3) the localised LSUBm scheme, in which all multi-spin correlations over all 
distinct locales on the lattice defined by m or fewer contiguous sites are retained. Another 
important feature of the method is that the bra and ket states are not always explicitly 
constrained to be Hermitian conjugates when we make such approximations, although the 
important Helmann-Feynman theorem is always preserved. We remark that the CCM pro- 
vides results in the infinite-lattice limit N — > oo from the outset. 

In this article, we wish to determine the spin-spin correlation functions |(s(0).s(r))| as 
a function of lattice separation \r\ for the (unfrustrated) square-lattice and (frustrated) 
triangular-lattice antiferromagnets. We must take into account the rotation of the local 
spin axes, although this proceeds in exactly the same manner as for the Hamiltonian above. 
Again, we remark that the manner in which high-order CCM may be solved has been 



discussed in Ref. [52|. The manner in which the ground-state CCM ket- and bra-state 
equations are solved is discussed in the Appendix to this article. In particular, the method 
by which "generalised expectation values" for a variety of one- and two-body spin operators 
may be obtained is explained in the appendices. 

III. RESULTS 

The results for the absolute values of the spin-spin correlation functions for the (un- 
frustrated) spin-half square-lattice antiferromagnet using the LSUBm approximation with 
m = {2, 4, 6, 8} are shown in Fig. HJ The LSUBm results are in very good mutual agreement 
for small lattice separations for different values of m. Indeed, LSUBm results are clearly 
converging with increasing levels of approximation. Furthermore, we see that LSUB6 and 
LSUB8 correspond reasonably well with each other up to r < 5. From Fig. [1] (as a "rule 
of thumb" only), we believe LSUBm results ought to provide reasonable results to lattice 
separations approximately of order m for the square lattice. 



CCM results are also compared to those results of exact diagonalisations (ED) [58| in Fig. 
CD However, we see that we obtain good correspondence only for very small lattice separations 
r < 3. (We note that N = 40 is at the current upper limit for ED of computational 
tractability.) However, the agreement between ED and CCM becomes rapidly worse for 
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FIG. 1: CCM results for the absolute values of the spin-spin correlation function for the spin-half 
square-lattice antiferromagnet using the LSUBm approximation with m = {2, 4,6, 8}. CCM results 



are compared to those results of exact diagonalisations (ED) with N = 36, 40 



58]. 



lattice separation r > 3. This 



disparity can be understood by considering results of quantum 



Monte Carlo (QMC) of Ref. [53] for the spin-spin correlation function, which were carried 
out for for much larger lattices of size of N = L x L, where L = {16,22,32}. Indeed, we 
find good correspondence qualitatively (i.e., compared by by eye only) between the LSUB8 
results of Fig. [1] and those QMC results for the largest lattice N = 32 x 32 of Fig. 3a in 
Ref. ||37|. Furthermore, the authors of Ref. noted a strong cusp-like behaviour was 
seen in the spin-spin correlation function at lattice separations of L/2 due to finite-lattice 
effects. This is also seen in Fig. CD for the ED results at r = 3, although this cusp may be 
exacerbated by a small "kink" in the "true" spin-spin correlation function that occurs at this 
point anyway (as seen in both CCM and QMC results). Interestingly, we see from Fig. 3b of 



Ref. [52] that the accuracy of QMC results for the spin-spin correlation functions becomes 
more problematic (again as a "rule of thumb" only) for separations of order r > L/2 for a 
given lattice size N = L x L. Hence, we might also reasonably expect that results of ED 
might only be good for separations of r < 3 for lattices of size N = 36 and n = 40. Indeed, 



we note that ED agrees wei 



1 with highly converged CCM results (and incidentally QMC 



results of Fig. 3a in Ref. 57j - again comparing by eye only) in this region. 
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FIG. 2: CCM results for the absolute values of the spin-spin correlation function for the spin-half 
triangular-lattice antiferromagnet using the LSUBm approximation with m = {2,3,4,5,6}. CCM 
results are compared to those results of exact diagonalisations (ED) 6^] with N = 36. 

The author of the QMC calculations in Ref. (5^] noted that: "to good approximation, 



the two-point function depends only on |r|." As also noted in Ref. 



571 ]. the long-range 



behaviour of correlation function gives a constant, thus indicating a long-range ordered 



ground state. We see this also at all levels of approximation as shown in Fig. [TJ Results 26] 
for the sublattice magnetisation have shown that extrapolation of the CCM LSUBm results 
gives good correspondence to the results of QMC, as noted above, where both methods 
indicate from the sublattice magnetisation that approximately 61% of the classical Neel 



order remains in the quantum limit. Finally, the author of Ref. [57j remarked that the 
approach to the constant for larger separations could be fitted by either an exponential 
or power-law form. We tried both for the LSUB8 data with r < 8 and found nothing to 
contradict this statement, although the power-law decay seemed to work better (i.e., had 
a lower residual error) that the exponential law. The residual sum-of-squares error for the 
power-law form was 0.0015 and for the exponential form was 0.0034. Fits of these forms to 



the LSUB8 data were carried out using the R statistics language: http://cran.r-project.org/ 
The results for the (frustrated) spin-half square-lattice antiferromagnet using the LSUBm 
approximation with m = {2, 3, 4, 5, 6} are shown in Fig. [2J Again, LSUBm results are clearly 
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converging with increasing levels of approximation, and again they are in very good mutual 
agreement for small lattice separations. CCM results are again compared to results of ED 
for A" = 36 (e.g., see Ref. 621]). and again good agreement is found for very small lattice 
separations (r < 3). Another strong cusp-like behaviour is seen in ED results for the spin- 
spin correlation function for r « 3. Again, it is highly likely that this is primarily due to 
finite-size effects. Indeed, we see that CCM results for r > 3 are quite different to those 
of ED. Furthermore, CCM LSUBm results initially decay strongly and then tend to a con- 
stant value for larger lattice separations. The behaviour of the correlation function for the 
triangular lattice is broadly similar to that of the square lattice. Again, extrapolation of 



LSUBm results |4l| for the sublattice magnetisation demonstrate that the spin-half trian- 
gular lattice antiferromagnet is Neel ordered. For LSUB6 with r < 5.5, we found that the 
initial decay of the correlation function for CCM data in this case could be fitted well by 
either an exponential or power-law form (for both forms, residual sum-of-squares: 0.0002). 



IV. CONCLUSIONS 

We have presented new formalism for the CCM in order to form "generalised" expectation 
values for a wide range of terms in the spin Hamiltonian. This is particularly useful for 



"high-order" CCM [63j. We have applied this approach in order to determine the spin-spin 
correlation functions for frustrated and unfrustrated 2D antiferromagnets, namely, the spin- 
half triangular- and square-lattice Heisenberg models. We found good correspondence of 
CCM with results of quantum Monte Carlo (QMC) for the square-lattice antiferromagnet. A 
strong cusp-like behaviour was noted in exact diagonalisation (ED) results for the correlation 
function at a lattice separation L/2 for a lattice of size N = L x L for both lattices. This is 
a finite-size effect. A cusp in the spin-spin correlation function was also seen in QMC results 
at r = 1/2 for a iV = L x L lattice. 

By contrast, this cusp was not seen in CCM results for either lattice at any value of r at 
any level of LSUBm approximation. Indeed, we note that CCM results are determined in 
the infinite-lattice limit (N — > oo) from the outset. We believe that (as a "rule of thumb") 
we obtain reasonable results for the correlation function up to separations of r « m for 
the square lattice - even despite the fact that we are using a purely "localised" LSUBm 
approximation scheme. LSUBm results for the triangular lattice are clearly converging, and 
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are in good mutual agreement for small lattice separations. For both lattices, CCM results 
for the correlation function decay strongly initially, although they become constant for larger 
values of r. Previous CCM results for the sublattice magnetisation 41] indicate that the 
spin-half Heisenberg models for the square and triangular lattices are Neel ordered, which 
agrees with the results of other approximate methods. 

It is still fair to say that QMC results still generally provide the most accurate results for 
the unfrustrated Heisenberg model on lattices of two spatial dimensions. However, QMC is 
severely restricted by the "sig n problem," which is a consequence of frustration at T = 0. 
The CCM has been shown 4l| to provides good results for the ground-state properties of the 
spin-half triangular lattice antiferromagnet such as the ground-state energy and sublattice 
magnetisation. Here we have shown that the CCM also provides reasonable results for the 
spin-spin correlation function of the spin-half triangular lattice antiferromagnet. Finally, 
we note that CCM results indicate that the behaviour of the correlation function for the 
triangular lattice is broadly similar to that of the square lattice. 
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APPENDIX A: CCM FORMALISM 



We begin the description of the application of the CCM by noting that it may be proven 
from Eqs. (J2J) and Q in a straightforward manner that the ket- and bra-state equations are 
thus given by 

($\Cj-e- s He s \$) = 0, V/ + ; (Al) 
(<£>\Se- s [H,C?]e s \<f>) = 0, V I ^ . (A2) 

The index I refers to a particular choice of cluster from the set of (Np) fundamental clusters 
that are distinct under the symmetries of the crystallographic lattice and the Hamiltonian 
and for a given approximation scheme at a given level of approximation. We note that these 
equations are equivalent to the minimization of the expectation value of H = (^>\H\^/) 
with respect to the CCM bra- and ket-state correlation coefficients {Si, Si}. We note that 
Eq. (|Aip is equivalent to 5H/5Si = 0, whereas Eq. (1A2|) is equivalent to 5H/5Si = 0. 
Furthermore, we note that Eq. flAip leads directly to simple form for the ground-state 
energy given by 

E 9 = E g ({Si}) = (*\e- s He 3 \*) . (A3) 

The full set {Si, Si} provides a complete description of the ground state. For instance, an 
arbitrary operator A will have a ground-state expectation value given as 

A = $\A\y) = (<5>\Se- s Ae s \<5>) = A ({Si, Si}) . (A4) 

The similarity transform of A is given by, 

A = e ~ s Ae s = A + [A, S] + h[A, S], S] + • • • . (A5) 



APPENDIX B: HIGH-ORDER CCM 



The manner is which the ground-state problem is solved to high-orders of approximation 



via computational is provided in Ref . [52J . However, we note that use new "high-order CCM 
operators" that are formed purely of spin-raising operators with respect to the model state. 
The model state is taken to be a state in which all spins point in the downwards ^-direction 
after some appropriate rotation of the local axes of the spins. This allows us to write one- 
and two-body spin operators in terms of these new high-order CCM operators. The problem 

11 



inherent in Eq. (1A1I) becomes one of pattern-matching the fundamental set of clusters in 
Cj to the terms in the similarity transformed version of the Hamiltonian e~ s He s \&) , where 
the 7-th such equation is given by 

Ei = {<!>\Cje- s He s \§} = ,V/^0 . (Bl) 

(Note that we assume that (<fr\CfCi~\<&) = 1 in the above equation). Specific terms in the 
Hamiltonian that may be used in the CCM code are: s z s z ; s + s z ; s~s z ; s z s + ; s z s~; s + s~; 
s~s + ; s + s + ; s~s~; s + ; s~; s z ; and, (s z ) 2 . We now "pattern-match" the C~ operators to 
those the relevant terms in the similarity transformed version of the Hamiltonian in order 
to form the CCM equations Ej — of Eq. fIBip at a given level of approximation. 

We now define the following new set of CCM bra-state correlation coefficients given by 
xj = Si and xj = Nb/N(1\)viSi and we assume again that (<&\CfCi~\<f>) = 1- Note that 
Nb is the number of Bravais lattice sites. Note also that for a given cluster I then vi 
is a symmetry factor which is dependent purely on the point-group symmetries (and not 
the translational symmetries) of the crystallographic lattice and that I is the number of 
spin operators. We note that the factors uj, N, Nb, and (Z!) never need to be explicitly 
determined. The CCM bra-state operator may thus be rewritten as 

Np 

S = 1 + NJ2 xiCj , (B2) 
i=i 

such that we have a particularly simple form for H, given by 

Np 

H = NJ2S: I E I , (B3) 

1=0 

where xq — 1. We note that the Eq is defined by Eq = jf(&\e~ s He s \&) (and, thus, Eq = 
jjE g ) and that E 1 is the J-th CCM ket-state equation defined by Eq. (jBTjl . The CCM ket- 
state equations are easily re-derived by taking the partial derivative of H/N with respect to 
xi, where 

«(,»)^„ (B4) 

OXi 

We now take the partial derivative of H/N with respect to Xi such that the bra-state 
equations take on a particularly simple form, given by 

M = ^ + g Ao) = £, . (B5) 

OX j OX i jz x OX i 
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APPENDIX C: GENERALIZED GROUND-STATE EXPECTATION VALUES 



The expectation value of a "generalized" spin operator that we shall call A may be treated 
in an analogous manner to that of the expectation value of the Hamiltonian, given by H . 
We write: 

A 7 = (§\Cje- s Ae s \<&) (CI) 

and with Cq = 1 . The similarity transform of A is defined by Eq. (1A4|) . We may treat a wide 
range of one- and two-body spin operations. However, unlike the Heisenberg Hamiltonian 
of Eq. ([1]), we do not constrain i and j in the two-body terms to be only nearest-neighbors 
or next-nearest-neighbors. For example, we consider here the spin-spin correlation func- 
tions. The expectation value of the generalized (spin) operator may again be written in a 
particularly simple form as: 

A = ($\A\y) 

N F 

= Nj^xji^lCye-SAe 3 ^} 
1=0 

N F 

^A = Nj^xjAj , (C2) 

7=0 

where Xo = 1 also and Aq = jj{Q\e~ s Ae s \&). The same code used to find ground-state 
equations may be used to find the generalized expectation values. Again we note that the 
index I in Eq. (1C2I) runs from zero to Np. Again, we note that factors such as Nb or uj etc. 
do not need to be determined explicitly because they cancel because of the definition of {xi} 



given above. The techniques need to achieve a computational solution have 



26 



3een_discussed 
3t 



32 



52j for 



extensively elsewhere, and the interested reader is referred to Refs. 
more information. However, clearly we see that the summation over all fundamental clusters 
involved in evaluating A is carried out readily once the ket- and bra-state equations of Eqs. 
( 1B4I) and (IB5I) have been solved for. 
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